Dielectric functions and collective excitations in MgB2 
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The frequency- and momentum-dependent dielectric function e(q, tj) as well as the energy loss 
function Im[-e~^(q, oj)] are calculated for intermetallic superconductor MgB2 by using two ah ini- 
tio methods: the plane-wave pseudopotential method and the tight-binding version of the LMTO 
method. We find two plasmon modes dispersing at energies ~ 2-8 eV and ~ 18-22 eV. The high 
energy plasmon results from a free electron like plasmon mode while the low energy collective exci- 
tation has its origin in a peculiar character of the band structure. Both plasmon modes demonstrate 
clearly anisotropic behaviour of both the peak position and the peak width. In particular, the low 
energy collective excitation has practically zero width in the direction perpendicular to boron layers 
and broadens in other directions. 
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After the discovery of superconductivity in the MgB2 
compound with the transition temperature Tc Z'dK 
101 much effort has been devoted to understanding the 
mechanism of the superconductivity [§ 11| as well as to 
studying different electronic and atomic characteristics 
of this compound. Among these characteristics are the 
superconducting gap p|JT2 Jl^ , the crystal structure and 
its influence on 0g,|l4jl6| , band structure ilj^Ml] 
and the Fermi surface U, lattice vibrations ||6|,|9|, 20,pi[ 
as well as thermodynamic and transport properties [ 22[ . 
Recently Voelker et al Q explored collective excitations 
very near the Fermi level by using a simple band struc- 
ture model and found the plasmon acoustic mode at very 
small momenta {q ~ 0.01 a.u.^^) and low energies {uj ~ 
0.01 eV). Here we study collective excitations in MgB2 
for different momenta and energies, namely, for q >0.1 
a.ur^ and w > 2 eV. We report first-principle calcula- 
tions for the real, ei(q, w), and imaginary, €2{(\,uj), part 
of the dielectric function as well as the energy loss func- 
tion Im[-e^^(q, cj)]. As a result of the calculation two 
plasmon modes and a few features arising from intcrband 
transitions are obtained. 

MgB2 crystallizes in the so-called AIB2 structure in 
which B atoms form graphitelike honeycomb layers that 
alternate with hexagonal layers of Mg atoms. The mag- 
nesium atoms are located at the center of hexagons 
formed by borons and donate their electrons to the boron 
planes. Similar to graphite MgB2 exhibits a strong 
anisotropy in the B-B lengths: the distance between the 
boron planes is significantly longer than in-plane B-B 
distance. We use this resemblance between graphite and 
MgB2 in order to clear up the origin of plasmon peaks in 
MgB2 by comparing the calculated energy loss function 
features with those obtained from EELS measurements 
for graphite and single-wall carbon nanotubes (SWCN) 
1^,^. To shed more light on the problem we also eval- 
uate the dielectric functions and the energy loss spectra 
for the MgB2 crystal structure with the Mg atoms re- 



moved (this hypothetical crystal structure is designated 
as 52). 

Information on the energy lost by electrons in their in- 
teractions with metals is carried by the dynamical struc- 
ture factor S{q,uj) which is related by the fluctuation- 
dissipation theorem to the energy loss function Im[- 
eQQ^(q, w)]. To calculate the inverse dielectric function 
we invoke the random phase approximation (RPA) where 
e^^(q, cj) is defined as (in symbolic form) = 1 + 
fc(l — x'^fc) x"; where Vc is the bare Coulomb poten- 
tial and x° is the density response function of the non- 
interacting electron system. The dielectric function is 
related to x'' as e = 1 — UcX*^- The energy loss function 
may be obtained by inverting the first matrix element of 
e that leads to neglecting short-range exchange and cor- 
relation effects or directly from when these effects are 
included. We have computed the energy loss function by 
using both of these approaches and found that the inclu- 
sion of the local field effects leads to negligible changes of 
both the width and energy of the plasmon peaks. In the 
calculation of the density response matrix Xqq' (q, i^) we 
have used two different first-principle methods: the plane 
wave pseudopotential method [E3| and the tight-binding 
version of the LMTO method [|2|,|||. 

In Figs, la and lb we show the evaluated band struc- 
ture of MgB2 and B2 along the symmetry directions. 
In general, these band structures are quite similar. The 
distinctions between them in the vicinity of the Fermi 
level (Ep) are due to the lower position of Ep relative 
to the a band in the TA direction for B2. The states of 
this band, which are of px.y symmetry, are degenerate in 
TA and their charge density is located in B layer show- 
ing a clear 2D character. This character leads to weaker 
interactions between the B layers and to smaller disper- 
sion of the a band along FA in B2 ■ The Pz band which 
is occupied at F in MgB2 becomes unoccupied at F in 
B2. In Figs. 2a and 2b we present the momentum de- 
pendence of the energy loss function in the F^, TK and 
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TM directions. In MgB2 we have found two plasmon 
modes. The higher collective excitation mode originates 
from the free electron like excitation mode with energy 
a;pi = 19.1 eV that corresponds to the electron density pa- 
rameter rs = l.82 a.u. of MgB2- This free electron like 
mode is transformed into two separated submodes 
and LOpi in a real crystal. One of them is very isotropic 
in the hexagonal plane (the TK and TM directions) and 
disperses linearly up for momenta q > 0.2 a.u., while an- 
other one, in the TA direction, has smaller energy and 
is nearly constant. Because of limitations of the calcula- 
tion methods and a large width of the energy loss peaks 
at small momenta we could not determine accurately the 
plasmon peak position in this region. Therefore we de- 
fine the plasmon energy at the T point by extrapola- 
tion of the computed plasmon dispersions at q> 0.2 a.u. 
This extrapolation results in 07^^=19.4 eV and Wpj=18.8 
eV in good agreement with the free electron gas value 
ujpi. The width, Api, of both these energy loss peaks 
decreases with the increasing momentum, and the A^i 
width deacreases faster than A^^. 

A different behavior is shown by the low-energy loss 
function peak which disperses linearly up in both the 
TA direction and in the hexagonal plane. In the TA 
direction, at q± ~ 0.1 a.u.~^, the peak is very narrow, 
^ 0.01 eV, that can be seen in the very small value of 
£2 in the energy interval around the peak position where 
£1=0 (Figs. 3a and 3b). In particular, for q±=0.12 a.u."^ 
this interval is between 1 and 5 eV (Fig. 3a) and the 
energy loss peak is located at 2.9 eV. On changing the 
momentum to the A point this interval becomes more 
narrow and moves to higher energies (Fig. 3c). At small 
momenta the first peak of Im[-£oo(q, w)] placed in the 
energy interval 0-1 eV is determined by intraband tran- 
sitions within the two a bands in the x, y plane around 
the TA direction, while the second peak located at 5.4 eV 
(Fig. 3a) is determined by the interband cr — tt transitions 
in the KM, AH, and AL directions. So, the low energy 
plasmon excitation corresponds to electron excitations in 
the a bands and one can define it as the a plasmon. In 
the hexagonal plane, in the TK direction, the low energy 
EELS peak broadens (Fig. 4a) and disperses linearly up 
on going from the F point to K. In the TM direction the 
plasmon peak disperses similar to that in the TK one 
showing a nearly ideal isotropy in the hexagonal plane. 
However, it becomes smaller and wider on going from F 
to M and vanishes finally at q ~ 0.8| TM \. Comparing 
the plasmon energies obtained from the LMTO and pseu- 
dopotcntial calculations one can find only a small differ- 
ence of ^ 0.1 eV between them. For instance, at q—0.2 
a.ur^ the LMTO ujp2 is slightly smaller than the pseu- 
dopotcntial one and vice versa for larger momenta. This 
slight difference results in different energy loss peak po- 
sitions at F: the extrapolation of both plasmon energies 
uj^2 ^p2 calculated for q> 0.1 a.u.^^ to the F point 



gives u;^2=l-8 eV, a;p|=2.0 eV (LMTO) and w^2=2-2 eV, 
uj^y=2A (pseudopotential). We estimate the accuracy of 
these values to be better than 0.2 eV. 

Besides two plasmon modes we have obtained four 
small features in Im[-£QQ^(q, lj)] that correspond to in- 
terband excitations. It is difhcult to find out what tran- 
sitions are responsible for these features, nevertheless we 
show them in Fig. 2a. One of them arises at q ~0.1a.M.~^ 
at an energy of w2 eV in the hexagonal plane, another 
one occurs at g > 0.4 a.u.~^ in the TM direction at an 
energy of ~ 4 eV and the other two small features arise 
in the TA direction for g=0.1-0.25 a.u.^^ at energies of 
10 eV and 13 eV, respectively. 

In Fig. 2b we show the momentum dependence of the 
energy loss function calculated for the hypothetical crys- 
tal structure B2. In general, the energy loss function in 
B2 shows relatively similar features to those in MgB2, 
though there are some important distinctions. In par- 
ticular, all collective excitations including two plasmon 
modes manifest a smaller dispersion in the F^ direction. 
This effect is a direct consequence of a weaker interac- 
tions between adjacent layers of boron in B2 compared 
to MgB2- Another distinction is that all features in the 
energy loss function in B2 are much clearer than those 
in MgB2 (Figs. 3a-3c and 4a-4c). One exception is the 
high energy plasmon mode. The third distinction is that 
B2 has more features in Im[-£gQ^(q, w)] than does MgB2- 
The low energy plasmon mode extrapolation to the F 
point gives ijjp2=^A eV which is ^ 2 eV higher than that 
in MgB2. This shift in energy is due to the higher en- 
ergy position of the second maximum of Im[£oo] (Fig. 3a) . 
While the position of the first peak of Im£ in B2 nearly 
coincides with that in MgB2 the second peak is moved 
by 2 eV to higher energies. Via the Hilbert transform 
(Kramers-Kronig relation) it also moves the node of £1 
to higher energy. Despite some quantitative distinctions 
between the energy loss functions in MgB2 and B2 one 
can conclude that mostly the features of the excitation 
spectrum of MgB2 can be derived, with the relevant cor- 
rections, from those of the hypothetical crystal i?2. 

The two plasmon modes similar to those obtained 
in MgB2 were also observed in EELS experiments 
for graphite and SWCN |26|j2^ ] which are even more 
anisotropic than MgB2. The upper plasmon mode which 
has larger energy in graphite and SWCN than in MgB2 
(near the F point w^f -21 eV (SWCN) and w^f -26 eV 
(graphite)) |^ also results from excitations of all valence 
electrons. The lower plasmon mode shows a linear 
dependence on momentum, like that in MgB2, with en- 
ergies c^pl ~5 eV (SWCN) and u^l ~6.5.eV (graphite) 
prf at F. But in contrast to SWCN and graphite where 
ujp2 represents the collective excitation of the 7r-electron 
system [p6|-p8[ in MgB2 this mode is a result of the col- 
lective excitation of the cr-electron system. The different 
origin of the low-energy plasmon mode in MgB2 / B2 and 
graphite/SWCN can be qualitatively understood from 
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the Fermi energy (Ep) position. In graphite the Fermi 
level is pinned by tt- electrons in the KH direction. In 
MgB2 and B2 with a smaller number of electrons per 
atom the Fermi level is pinned by cr-electrons. So, low- 
energy excitations in MgB2 and in B2 are expected to 
be derived from the cr-band electrons. 

In conclusion, we have performed the first-principle 
calculations of the dielectric functions ei(q, w) and 
e2(q, w) as well as the energy loss function Im[-e~^(q, oj)]. 
The calculations reveal the two plasmon modes in MgB2 
and B2 and a few interband collective excitations. The 
low energy plasmon mode corresponding to the excita- 
tions of electrons in the a bands shows a very anisotropic 
behavior of the peak width. The energy loss spectrum of 
MgB2 can be derived, with the relevant corrections, from 
that of the hypothetical crystal structure B2 ■ 
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After the submission of this paper first-principles cal- 
culations of collective excitations in MgB2 have also been 
presented by Wei Ku et al. |29] . 
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FIG. 1. Calculated energy band structure of a) MgB2 and 
b) B2 along the symmetry directions, a and n denote the 
boron bands oi px,y and pz character, respectively. 



FIG. 2. (color) Dispersion of the plasmon modes (filled cir- 
cles) and of the features (triangles) for a) MgB2 and b) B2 
along the symmetry directions: TA (green), TK (red), and 
FM (blue). 
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FIG. 3. (color) Real (dashed lines) and imaginary (solid 
lines) parts of the dielectric function eGG{q,u>) along the TA 
direction for G=0 and q=(0,0,Q)(27r/c): a) Q=l/8, b) q=1/4, 
c) Q=l/2 (A point). The energy loss function is shown in the 
insert. The red (blue) lines represent MgB2 (B2). 



FIG. 4. (color) Real (dashed lines) and imaginary (solid 
lines) parts of the dielectric function eGG(q, oj) along the TK 
direction for G=0 and q=(Q,0,0)(27r/a): a) a=l/9, b) Q=l/3, 
c)a=2/3 (K point). The energy loss function is shown in the 
insert. The red (blue) lines represent MgB2 (B2). 
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Fig. 1b 
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